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Abstract 

We study the linear resistance at the Kosterlitz-Thouless transition by Monte 
Carlo simulation of vortex dynamics. Finite size scaling analysis of our data 
show excellent agreement with scaling properties of the Kosterlitz-Thouless 
transition. We also compare our results for the linear resistance with experi- 
ments. By adjusting the vortex chemical potential to an optimum value, the 
resistance at temperatures above the transition temperature agrees well with 
experiments over many decades. 

PACS numbers: 05.70Jk, 64.60Cn, 74.20De 
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Scaling of current-voltage characteristics is useful both in theoretical and experimental 
analysis of superconductors. One useful way to calculate linear and nonlinear resistance 
accurately is by computer simulation, and scaling analysis is crucial to extract the critical 
properties. The linear resistance and nonlinear current- volt age characteristic are key quanti- 
ties in the recently much studied vortex glass state in disordered superconductors ][]]. At the 
vortex glass transition the nonlinear current- volt age characteristic is universal, and the lin- 
ear resistance vanishes, signaling true superconductivity. The current-voltage characteristic 
is also a key quantity of the zero-field superconducting transition in two dimensions, which 
is a Kosterlitz-Thouless (KT) transition . The resistance curve above the KT transition 
is given by an (unknown) universal function of a reduced Coulomb gas temperature variable, 
allowing experimental data for different 2D systems to collapse on the same curve, as shown 
by Minnhagen [f|]. Remarkably, also certain 3D high temperature superconductors, with 
weakly coupled layers, seem to fit the same curve, except very close to T c [||J)|. The term 
"universality" for the resistance curve is used in the same sense as e.g. the universality of 
the BCS gap function, rather than in the sense of critical phenomena. The form of the 
resistance curve has been accurately constructed from such experiments, but an accurate 
calculation is lacking. KT theory gives an approximate expression for the resistance, valid 
in the KT critical region, but which cannot accurately fit the experimental curve outside the 
KT critical region (this will be seen below). One calculation going beyond the approximate 
KT resistance expression was a simulation of the sine-Gordon formulation of the Coulomb 
gas |7j]. An indirect connection was found to experiments, but a direct calculation of the 
resistance curve is lacking. 

In this paper we report results of Monte Carlo (MC) simulation and finite size scaling 
analysis of the linear resistance, from equilibrium vortex dynamics of the 2D Coulomb gas 
model for the KT transition. We first use finite size scaling analysis to extract critical 
properties from our MC data for the linear resistance, and we find excellent agreement with 
scaling properties of the KT transition. Secondly, we show that under certain conditions, our 
results for the linear resistance can closely fit the experimental curve above the transition 
temperature. 

Let us first consider scaling properties of the linear resistance at a superconducting 
transition driven by vortex fluctuations, like the KT transition or the vortex glass transition. 
At the critical temperature T c the correlation length diverges: £ — > oo. Dynamic quantities 
like the resistance depend on the correlation time r, which diverges as r ~ f , where z is 
the dynamical exponent. Voltage is related to the rate of phase slip of the superconducting 
order parameter, V = h/(2e)<p ~ £,~ z - The linear resistance is related to the equilibrium 
voltage fluctuations according to the Nyquist formula H-H] 



1 r+oo 

R =^fl dt (V(t)V(0)), (1) 



oc 



and thus the resistance scales as R ~ £~ z . 

Now we specialize to two dimensions. The superconducting transition is now a KT 
transition, where confined neutral vortex pairs start to deconfine, so that "free" vortices 



appear. Here we will assume the value z = 2 pi|J12|| , which describes free vortex diffusion. 
Hence the resistance is expected to scale like 

R-r 2 . (2) 
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This scaling form, here obtained from general scaling arguments, is the same as from the 
usual KT theory (see below). Simulations are done in finite systems, and the diverging 
correlation length at T c is cut off at the system size L. Precisely at T c the resistance is 
therefore expected to obey R ~ L~ 2 , or L 2 R = A, where A is a (nonuniversal) constant. 

As an aside comment, we note that another way of introducing a finite length scale at 
the transition is to apply a finite current density J. The finite current length £j, defined 
by £jJ = ksT is the length beyond which the diverging correlations at the transition 
are driven out of equilibrium by the Lorentz force. The resulting nonlinear response obeys 
E = R(J)J ~ ij z J ~ J z+1 = J 3 , which is the usual KT result, here obtained within a 
general scaling approach. 

We can extract the leading finite size correction to scaling of the resistance by using 
known results from KT theory. The correlation length £ diverges as H 



f- 2 



g 



(3) 



eT 

for T — > T c + , where g is a constant, 1/e is a dielectric function proportional to the renor 



malized superfluid density, and z is the fugacity. The Bardeen-Stephen relation [JT3 
R/Rn = 27i( 2 riF, ties the resistance to the density of free vortices rip, where Rn is the 
normal state resistance and ( the Ginsburg-Landau length, and np is in turn related to £ 
according to the Debye-Htickel relation [|J, = 2-RnpjeT . Eliminating rip gives back the 
scaling relation, R ~ £ -2 , in agreement with Eq. (El). The leading correction to scaling is 



known for 1/e JT4|: l/e(L) = l/e(oo)(l + l/(21n(L) + C)). Combining these relations gives 
again a finite size scaling relation for R, similar to L 2 R = A, but now containing the leading 
logarithmic correction: 

where C is an unknown constant. We also need the standard KT form of the resistance [HI. 



Using the temperature dependence of e from KT theory, £ becomes £ ~ exp( const j\/T — T c 
which gives 



R ~ g-const/v/T^ (5) 

expected to be valid in the KT critical scaling region. 

In the simulation we use the 2D Coulomb gas model. We consider a square lattice 
with N = L 2 sites with periodic boundary conditions. The grand partition function is 
Z = Tr{ n .yexp(-H/T), with 

H= \Y, n i n J G ij ~ V n l ( 6 ) 

i,j j 

where the trace is over rii = 0, ±1, where rii = means no vortex and = ±1 means one 
vortex of vorticity ± at site i. The parameter ix = —E c is a "vortex chemical potential" 
and E c is the vortex core energy. We only allow configurations of zero total vorticity. The 
vortex-vortex interaction is a lattice Green's function, G(r) = (2tt/N) ^k^o ex P(^k- r)/(4 — 
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2cos& £C — 2 008^), which varies like — lnr on large distances (the lattice spacing is set to 

0=1). 

Now we will describe the simulation. Following Ref. | 12"f , our MC algorithm consists of 



attempts to insert near neighbor (n = +l,n = —1) pairs, on randomly chosen lattice posi- 
tions, and with random orientation, and we test for acceptance according to the Metropolis 
algorithm. We typically discard 10 4 initial sweeps through the system to approach equilib- 
rium, and then take measurements during up to 10 6 sweeps close to T c , and fewer away from 
T c . We compute the linear resistance in the simulation by evaluating the Nyquist formula, 
Eq. ([I]), as a sum over discrete time steps PJTH[|. The voltage V{t) is given by ±1 (in suitable 



units), if at MC time t insertion of a vortex pair with orientation H — ( — h) is accepted, and 
V (t) = otherwise. To equate MC time to real time assumes heavily overdamped dynamics, 
which should be fulfilled near the transition. 

Now we turn to the results. In Figure |l] we determine the transition temperature from 
finite size scaling of MC data for the linear resistance. The vortex chemical potential was 
here chosen to p = 0, which gives the Coulomb gas corresponding to the 2D XY model in 
the Villain approximation. According to Eq. (|2]) data for L 2 R for different lattice sizes L 
should be independent of system size at T = T c . This means that T c is where the curves in 
the figure cross. Note that this scaling procedure assumes initial knowledge of the dynamical 
exponent, z — 2. Fig. |l] (a) is based on Eq. (|2|) and does not contain corrections to scaling. 
Fig. [l| (b) is based on Eq. (f|) which contains the lowest logarithmic correction to scaling. 
The constant C in Eq. (|j) was adjusted to as closely as possible have all curves cross at 
a single point. Comparing (a) and (b) shows that including the correction significantly 
improves the fit, and gives the estimate T c « 0.218 for p = 0. Error bars on the data points 
have been left out because they are much smaller than the symbol size. 

Figures |T] (a) and (b) show that the resistance obeys the expected scaling form R ~ L~ 2 
at the transition, T = T c . This also means that in the superconducting phase below T c , the 
resistance scales to zero with system size. This is possible in the Nyquist formula (|I]), despite 
that the equal-time voltage correlation, (V 2 (0)), gives a finite, roughly size independent, 
contribution to the resistance at all finite temperatures, also below T c . What happens is 
that below T c the unequal-time correlations on average cancel the equal-time correlations. 
This cancellation is incomplete above T c , and in finite systems also below. Within the KT 
vortex unbinding picture, below T c all vortices are bound in neutral pairs. Creation of a 
pair gives a burst of voltage that contributes to the equal-time voltage correlation, but if 
the pair is later annihilated a canceling unequal-time correlation results. Above T c some 
pairs unbind, the cancellation becomes incomplete, and the resistance is finite. 

We wish to compare the present approach of using the linear resistance to locate the 
critical point, to the often used method of locating the universal jump in the superfluid 
density p s . In the Coulomb gas p s corresponds to the dielectric response function p s ~ 
l/e(A; — > 0), where l/e(k) = 1 — 27r/(/c 2 TiV)(|?7.(k)| 2 ), where n(k) the the Fourier transform 
of the vortex density. The strong size dependence of R ~ L 2 produces much stronger splay 
of data away from T c , than the universal jump which has size dependence p s ~ L°. This 
makes locating T c easier. Furthermore, use of l/e(k) is in practice limited to the smallest 
available nonzero wave vector, k = 2tt/L, while R is for k = 0. 

Now we will compare resistance data to the experimental curve over a finite temperature 
interval above T c . In order to do this we first employ a crucial empirical correction to our 
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MC results. Our discrete MC dynamics fails in the limit of high temperature, because the 
voltage fluctuations saturate at high T when most MC moves are accepted, which incorrectly 
makes the resistance vanish like R ~ 1/T — > according to Eq. ([I]). Instead of approaching 
a constant at high T, we assume an additional factor of T in the vortex velocity squared, as 
for the temperature dependence of a diffusion process according to the Einstein relation || . 
This correction cancels the factor 1/T in the Nyquist formula, and makes R approach a 
constant at high T. The correction is negligible in the critical region, compared to the 
strong exponential temperature dependence in Eq. (|^). 

Figure ^ (a) shows MC data for the resistance curve, together with the experimental 
curve. Note that the temperature scale is here the Coulomb gas temperature, T = T , 
which differs from the physical temperature by the temperature dependence of the unrenor- 
malized superfTuid density p s o: 1/T CG = 27rp< j o/^B^ physlcal 0- The only parameters involved 
in the plot are T c , which is determined separately by the above scaling approach, and the 
normalization constant Rn, which is adjusted to obtain agreement with the experimental 
curve at high T. MC data for fi = is seen to deviate strongly from the experimental curve. 
To get better agreement we tried different values of the vortex chemical potential fi. We 
obtain best agreement for the value /i = —0.8. The critical temperature for /i = —0.8 is 
very close to the upper limit T c = 1/4, and the value T c = 1/4 is used in the plot. Further 
decreasing \i gives a worse fit. Figure | (b) is a KT plot with 1/(T/T C - l)" 1 / 2 on the x-axis. 
The KT form given by Eq. ([5]) corresponds to a straight line in this plot. 

Let us discuss this result. By adjusting the vortex chemical potential to /i = —0.8, our 
data for the resistance is in close agreement with the experimental curve over many decades, 
and over a broad temperature range. The agreement is not perfect, as indicated by the finite 
size effects near T c in Fig. ||, but the agreement is clearly better than a fit to KT theory 
outside the KT critical region. The value \i = —0.8 seems perhaps unexpectedly small, 
compared to e.g. the 2D XY model which has fi ~ 0. The corresponding Coulomb gas model 
has very small vortex density, and it is very close to the KT line of fixed points. One possible 
interpretation is that this is actually an effective vortex model with renormalized parameters. 
For example, vortex crystallization in superconducting films was suggested recently by Gabay 
and Kapitulnik (jTjJ. Our Coulomb gas model could describe logarithmically interacting 
vacancies with big core energy. 

We have also studied the effect of adding random vortex pinning, modeled as a random 
site energy for the vortices. This shifts T c downwards, but gives exactly the same scaling 
properties at the new T c . This directly verifies the expectation that weak disorder is indeed 
irrelevant at the KT transition, in agreement with the Harris criterion, v > 2/d flT6| , |T7 . 



In summary, our results from finite size scaling of MC data for the linear resistance of 
the 2D Coulomb gas are in excellent agreement with scaling properties of the KT transition. 
This shows that Monte Carlo dynamics is in the right dynamical universality class, and 
allows accurate determination of the critical temperature. Furthermore, by tuning the vortex 
chemical potential, our results can closely reproduce the experimental universal resistance 
curve over many decades, and the agreement outside the KT critical region is much better 
than a simple fit to the KT theory. 

We acknowledge stimulating discussions with S. M. Girvin, P. Minnhagen, H. J. Jensen, 
and S. Teitel. This work was supported by grants from the Swedish Natural Science Research 
Council (NFR). 
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FIGURES 



FIG. 1. Location of the critical temperature by finite size scaling of Monte Carlo data for the 
linear resistance R. The critical temperature is located where data for different system sizes L 
intersect. In (a) is shown data without the logarithmic correction to scaling, and in (b) is the same 
data including the logarithmic correction. 

FIG. 2. Resistance curves from Monte Carlo simulations for different values of the vortex 
chemical potential fi, together with the experimental curve. The resistance is plotted vs. T in (a), 
and vs. (T/T c — l) -1 / 2 in (b) as suggested by Eq. (||). 
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